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Abstract 

We calculate the CP averaged branching ratios and CP-violating asymmetries for rjr], rjrj' 

and 'q'r]' decays by employing the perturbative QCD (pQCD) factorization approach. The 
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for the considered three decay modes, which can be tested by the future B meson experiments. 
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I. INTRODUCTION 



The two-body charmless B meson decays provide a good place for testing the stan- 
dard model (SM), studying CP violation of B meson system, exploring the rich quantum 
chromodynamics (QCD) of strong interaction and searching for the signal or evidence of 
new physics beyond the SM Up to now, many B — »• M1M2 decays (where Mj refers 
to the light pseudo-scalar or vector mesons ) have been measured experimentally [2I, 0] 



with good accuracy, ca^ 
tion (QCDF) approach 



as 



10,1111,1121 

Among various B 
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culated and studied phenomenologically in the QCD factoriza- 
1, 1, 11 , B| , the perturbative QCD (pQCD) factorization approach 
15 . l^, or the soft collinear effective theory (SCET) (l7l |. 
M1M2 decay channels, the decays involving the isosinglet rj or rj' 
mesons in the final state have been studied extensively during the past decade because 
of the so-called Krj' puzzle or other special features. At present, we still do not know 
how large is the gluonic content of the rj' meson, and how to calculate reliably the gluonic 
contributions to the decay modes involving t]' meson. It is still difficult to explain the 
observed pattern of the branching ratios for B — > K^*^r]^'^ decays. 



In pQCD factorization approach, the B 



pr]'^'^ and ni]^'^ decays have been 



studied in Refs. [13|, UM llS 



In this paper, we would like to calculate the branching 
ratios and CP asymmetries for the three B rjrj, r]r]' and r]'r]' decays by employing the 



21I . Besides the 



low energy effective Hamiltonian [18| and the pQCD approach |19l . |2C 
usual factorizable contributions, we here are able to evaluate the non-factorizable and the 
annihilation contributions to these decays. On the experimental side, the CP-averaged 
branching ratios of i? — *■ rj^'^rj'^'^ decays have been measured very recently [22| in units of 
10-6 (upper limits at 90% C.L.): 

(1) 

(2) 
(3) 

The accuracy of the data is still low, only the upper limits can be used to compare with 
the theoretical predictions. 

This paper is organized as follows. In Sec. [TTl we give a brief review for the pQCD fac- 
torization approach. In Sec. IIIIl we calculate analytically the related Feynman diagrams 
and present the various decay amplitudes for the studied decay modes. In Sec. [TV] we 
show the numerical results for the branching ratios and CP asymmetries of i? ^ rj^'^^'^ 
decays and compare them with the measured values or the theoretical predictions in 
QCDF approach. The summary and some discussions are included in the final section. 
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II. THEORETICAL FRAMEWORK 



For the non-leptonic B decays, the dominant theoretical uncertainty comes from the 
evaluation of the hadronic matrix element < MiM2\0i\B >. Now there are two popular 
factorization approaches being used to calculate the matrix elements: the QCDF approach 
[3] and the pQCD approach [l^, 20, 21|. The pQCD approach has been developed earlier 
from the QCD hard-scattering approach [2l|. Some elements of this approach are also 
present in the QCD factorization approach^ll]. The two major differences between these 
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two approaches are (a) the form factors are calculable perturbatively in pQCD approach, 
but taken as the input parameters extracted from other experimental measurements in 
the QCDF approach; and (b) the annihilation contributions are calculable and play an 
important role in producing CP violation for the considered decay modes in pQCD ap- 
proach, but it could not be evaluated reliably in QCDF approach. Of course, one should 
remember that the assumptions behind the pQCD approach, specifically the possibility 
to calculate the form factors perturbatively, are still under discussion [23|. More efforts 
are needed to clarify these problems. 

In pQCD approach, the decay amphtude is separated into soft, hard, and harder dy- 
namics characterized by different energy scales {t,mh, Myy). It is conceptually written as 
the convolution, 

AiB^M^M2) ~ Jd%d^k2d% Ti[C{t)<l>Bih)<l>M,{k2)<^M,{k3)H{h,k2,h,t)], (4) 

where fc^'s are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C{t) is the Wilson coefficient which results from the radiative 
corrections at short distance. In the above convolution, C{t) includes the harder dynamics 
at larger scale than Mb scale and describes the evolution of local 4-Fermi operators from 
rriw (the W boson mass) down to t ~ 0{\/AMb) scale, where A = MB—rrib. The function 
H{ki, /c2, fcs, t) is the hard part and can be calculated perturbatively. The function $m is 
the wave function which describes hadronization of the quark and anti-quark to the meson 
M. While the function H depends on the process considered, the wave function $jv/ is 
independent of the specific process. Using the wave functions determined from other well 
measured processes, one can make quantitative predictions here. 

Since the b quark is rather heavy we consider the B meson at rest for simplicity. It 
is convenient to use light-cone coordinate (p^,P~,Pt) to describe the meson's momenta: 
p± = (p° ±p^)/^/2 and pr = 

Using the light-cone coordinates the B meson and the two final state meson momenta 
can be written as 

Mr Mr Mr 

Pi = -|(l,l,Or), P2 = ^(1,0,Ot), P3 = ^(0,1,Ot), (5) 

respectively, here the light meson masses have been neglected. Putting the light (anti-) 
quark momenta in B, rj' and 1] mesons as ki, k2, and k^, respectively, we can choose 

fci = (xiP+,0,kiT), A;2 = (a;2P2+,0,k2T), fcg = (0, X3P3", kar). (6) 

Then, for B rjr]' decay for example, the integration over fcf, fc^, and k^ in eq.(jl]) will 
lead to 

MB ^ ,„/) ~ 

■Tl [C(()4.B(ii,6,)t(r,')(i2,62)4„(i3,63)i/(i„*.,()Si(i.)e"""] . (7) 

where bi is the conjugate space coordinate of kix, and t is the largest energy scale in 
function H{xi, hi, t). The large logarithms ln{mw/t) are included in the Wilson coefficients 
C{t). The large double logarithms (In^Xj) on the longitudinal direction are summed by 
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the threshold resummation [24'], and they lead to St{xi) which smears the end-point 
singularities on Xj. The last term, e""^^*-*, is the Sudakov form factor which suppresses the 
soft dynamics effectively [1^. Thus it makes the perturbative calculation of the hard 
part H applicable at intermediate scale, i.e., Mb scale. We will calculate analytically the 
function H{xi,bi,t) for the considered decays in the first order in as expansion and give 
the convoluted amplitudes in next section. 



Wilson Coefficients 



For the two-body charmless B meson decays, the related weak effective Hamiltonian 
Heff can be written as [181 



V2 



10 



i=3 



(8) 



where Ci{^) are Wilson coefficients at the renormalization scale /x and Oi are the four- 
fermion operators for the case oih d transition. 



O9 



daJ^LUfS ■ UfS'Jf.Lba , 



Ob 



daYLbp ■ Eg' Q'^li^Lq'a > 
\daYLbp ■ Eg' eq'q'p-ff,Rq'a 



(9) 



^daYLba ■ Eg' eq'q'f^^t^Lq'/s > Olo = ^d^-f'^Lb^ ■ Eg' ^q' q'f^-f f,Lq'^ 



where a and P are the SU{3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = {1 + 75). The sum over q' runs over the 
quark fields that are active at the scale fi = 0{mb), i.e., q'e{u,d, s,c,b}. For the Wilson 
coefficients Ci(/i) (i = 1, . . . , 10), we will use the leading order (LO) expressions, although 
the next-to-leading order (NLO) results already exist in the literature [l8|. This is the 
consistent way to cancel the explicit /x dependence in the theoretical formulae. For the 
renormalization group evolution of the Wilson coefficients from higher scale to lower scale, 
we use the formulae as given in Ref.Q directly. 



B. Wave Functions 



In the resummation procedures, the B meson is treated as a heavy-light system. In 
general, the B meson light-cone matrix element can be decomposed as |j, l26[ 



' d^z 

w? 

i 



iki z 



{OMO)d^{z)\B{pB)) 



'>B 



+ "^b)75 



i-ij 

7r< 



/3o 



(10) 



where n = (1,0, Ot), and v = (0, 1, Ot) are the unit vectors pointing to the plus and 
minus directions, respectively. From the above equation, one can see that there are 
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two Lorentz structures in the B meson distribution amplitude (DA). They obey to the 
following normalization conditions 



(27r)4^"' 2v^' J (27r 

In general, one should consider these two Lorentz structures in calculations of B meson 
decays. However, it can be argued that the contribution of (ps is numerically small 



27|, |28||, thus its contribution can be numerically neglected. Therefore, we only consider 



the contribution of Lorentz structure 

'^B = ^=(^B + "^B)750B(kl), (12) 
V ^^Vc 

in our calculation. We use the same wave functions as in Refs. [sl, 0, [28| . Through out 
this paper, we use the light-cone coordinates to write the four momentum as {k^ , , k^). 
In the next section, we will see that the hard part is always independent of one of the 
k^ and/or fcf , if we make some approximations. The B meson wave function is then the 
function of variable k^ (or kf) and k^. 

MK, ki) = ! dkt<P{ki, k^, ki). (13) 



The wave function for dd components of r/'-'^ meson are given as 14 



%JP, X, C) = Mi,{x) + mr>L(^) + - V ■ nW^Jx) , (14) 



where P and x are the momentum and the momentum fraction of r/^j respectively, while 



'Vdd' ^Vdd ^Vdd ^^present the axial vector, pseudoscalar and tensor components of 



the wave function respectively. Following Ref. [ij], we here also assume that the wave 
function of t^^j is same as the vr wave function based on SU(3) flavor symmetry. The 
parameter ( is either +1 or —1 depending on the assignment of the momentum fraction 

X. 

The transverse momentum fc-*" is usually converted to the b parameter by Fourier 
transformation. The initial conditions of the function (j)i{x) with i = [B,-)],-)]') are of 
non-perturbative origin, satisfying the normalization 

.1 ^ 

0,(x, b = 0)dx = —j=fi , (15) 
2VD 

with fi the meson decay constants. 



III. PERTURBATIVE CALCULATIONS 



In this section, we will calculate and show the decay amplitude for each diagram 
including wave functions. The hard part H{t) involves the four quark operators and the 
necessary hard gluon connecting the four quark operator and the spectator quark. We first 
consider B rjr]^'^ decay mode as an example, and then extend our study to 5 — >■ ?]'?]' 
decay. 
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A. Decay amplitudes 



Similar to the B txt]^'^ decays in 16j , there are 8 type diagrams contributing to 
the B — > r/r/*^') decays, as illustrated in Figure 1. We first calculate the usual factorizable 
diagrams (a) and (b). Operators Oi, O2, O3, O4, O9, and Oio are {V — A){y — A) currents, 
the sum of their amplitudes is given as 

/•i 1*00 

Feri = STcCF'm'^ / dxidxs / 61^6163^63 61) 

Jo Jo 

X { [(1 + xs)4>^{xs, 63) + (1 - 2x3)r^{(pP{x3, 63) + ^lixs, 63))] 

■as{tl) K{xi, X3, 61, 63) exp[-Sabitl)] 
+2rn4)^{x3,b3)as(tl)he{x3, xi,b3,bi) exp[-Sab(tl)]} . (16) 

where = mo/m^; Cp = 4/3 is a color factor. The function h^, the scales and the 
Sudakov factors Sab are displayed in Appendix [XI 

The form factors of B to rj decay, Fq^^{0), can thus be extracted from the expression 
in Eq. ([16]), that is 

Fo'^^'iq' = 0) = Ff-''(g' = 0) = F,,/m|, (17) 



which is identical with that defined in Ref. [2j 

The operators O5, Oq, O7, and Os have a structure of (V — A) (V + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the axial- vector part of (V + A) current contribute to the pseudo-sealer meson 
production, {ri\V - A\B){ri'\V + A\0) = -{r]\V - A\B){ri'\V - A\0), that is 

F^ri^ = -Ferj ■ (18) 

For other cases, we need to do Fierz transformation for the corresponding operators to 
get right fiavor and color structure for factorization to work. We may get {S — P){S + P) 
operators from {V — A){V + A) ones. For these {S — P){S + P) operators. Fig. 1(a) and 
1(b) give 



F^^ = 167rCi7m^r (,) / dxidxs I 61^6163^63 0b(xi, 6 







X { [0^(X3, 63) + r^((2 + X3)0^(X3, 63) - X30* (X3, 63))] 

■as{tl)K{xi, X3, 61, 63) exp[-Sabitl)] 
+ [xi0,,(x3,63) - 2(xi - l)r^0;^(x3,63)] 

■as{tl)K{x3,Xi,b3,bi)exp[-Sabitl)]} . (19) 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using S function 5(63 — 61), leaving only 
integration of 61 and 62. For the {V — A){V — A) operators, the result is 

16v^6 

Me-q = — - — nCpm^ / dxidx2dx3 / 61(^6162^62 0s (xi, 61) 0^(,) (x2, 62) 
<J Jo Jo 

X { [2x3r^0* (x3, 61) - X30^(a;3, 61)] 

■as{tf)hf{xi,X2, X3, 61, 62) exp[-Scditf)]} . (20) 
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V B 




(b) 
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V B 



(d) 
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B 




B 




FIG. 1: Typical Feynman diagrams contributing to the B ijrj' decays, where diagram (a) 

Q 

and (b) contribute to the B ^ rj form factor Fq^ . 

For the {V — A){V + A) operators the formulae are different. Here we have two kinds 
of contributions from {V -A)(y + A) operators: M^^^ and M[^^ is for the {V -A)(y + A) 
and (5 — P){S + P) type operators respectively: 



M^^ = 0, M^^ = -M 

eri en < 



er) 



'■er) ■ 



(21) 



The factorizable annihilation diagrams 1(g) and 1(h) involve only rj and rj^'^ wave 
functions. There are also three kinds of decay amplitudes for these two diagrams. Far, is 
for {V — A){V — A) type operators, is for {V — A){V + A) type operators, while 
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is for {S — P){S + P) type operators. 

^arl'^^ar, = -ST^Cpm^ I dX2 dx^ / 62C?&2&3C^&3 { [a;30^(x3, fegj^^w (^2, 62) 

Jo Jo 

+2r^r^(o{{x3 + l)(/)J(a;3, 63) + (x-^ - 1)0* (0:3, 63))0^(,) (x2, 62) 

■as(tl)ha{x2, X3, 62, ^3) exp[-5'g/i(tg)] 

- [a:20,,(x3,63)0,,(,)(a;2,&2) 

+%V((^2 + 1)0^W(X2, 62) + (X2 - l)0^(,)(a;2, b2))(l)^{x3, 63) 

•a.(i^)/i„(x3, X2, 63, &2) exp[-5^^(i^)]} (22) 



^P2 
at] 



pi poo 

— 167rCi?m5 / dx2dx3 / b2db2b3db3 
Jo Jo 

X I [a:3'^r;(0j'(a^3, ^s) " 0^(2^3, &3))0,,(') (2^2, ^2) + 2r^(,)0^(x3, 63)^50(^2, ^2)] 



■as{tl)ha{x2, X3, 62, ^3) exp[-S'g/i(ig)] 
+ [2^2^) (^J/) (a:;2, ^2) - 0*(o (2^2, b2))(f)r,{x3, bs) + 2r^0^(,) (x2, b2)(f)^{x3, 63)] 
■as{ti)ha{x3, X2, &3, ^2) exp[-S'g/j(tg)]} . 



(23) 



For the non-factorizablc annihilation diagrams 1(e) and 1(f), again all three wave 
functions are involved. Here we have three kinds of contributions. M^^, M^^ and M^^ 
describe the contributions from {V - A){V - A), {V - A){V + A) and {S - P)iS + P) 
type operators respectively: 



M. 



arj 



— ^^^7rCFm% f dx\dx2dx3 f bidbib2db2(f)B{xi,bi) 
3 Jo Jo 

X { - {X2(l)^{x3,b2)(t),,{0{x2,b2) 

+rr,r^(o {x2 + X3 + 2)0j,) {x2, 62) + {^2 - a;3)0^(/) ix2, 62) 0^(2:3, ^2) 

+r^r^(/) {x2 - X3)0j,)(a;3, 62) + (2^2 + 2^3 - 2)0*(,)(a:;2, ^2) ^^a^s, 62)} 
•Q;s(tJ)/iJ(a;i, X2, X3, bi, 62) exp[-5'e/(i^)] 

+ {X3(f)r,ix3, &2)0^(') (3:^2, ^2) 



{x2 + a;3)05,) (a;2, 62) + (2:3 - a;2)0* w (2:2, ^2) <P^{x3, ^2) 



+r^r^(,) (xs - a;2)05,) (a;2, 62) + {x2 + 2^3)0* w (2^2, ^2) (plixs, ^2) } 
•Q;s(t))/iJ(a:;i,X2,X3,6i,62)exp[-5'e/(i))]} , (24) 
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M, 



PI 



-^^^nCpm'^ f dxidx2dx3 f bidbib2db2(pB{xi,bi) 
3 Jo Jo 

X { [{X3 - 2)n^0^(,)(x2,&2)(0j'(a;3,&2) + 0*(X3,62)) - - 2)r^(,) 0^ (xs , 62 

(x2, 62) + (a;2, ^2)) ■ as{t^f)h^f ixi,X2, xg, 61, 62) exp[-S'e/(tJ)] 
- [a;3r^0^(,)(x2,62)(0j'(a;3,&2) + 0*^(x3,62)) 
-X2V)0r,(a;3,&2)(</>5o (3^2,^2) + (a;2,&2)) 
a;2, 3:3, ^i, ^2) exp[-S'e/(t))] } , 



(25) 



M, 



P2 

arj 



-nCprrL^B / dxidx2dxz / bidbib2db2(f)Bixi,bi) 
Jo Jo 



X {{a;30r,(a;3,&2)0,,(')(a;2,&2) 



+r„r, 



{X2 + X3 + 2)0^(,) (X2, 62) + (3:3 - X2)(j)li^,) ix2, 62) (I)^ix3, 62' 

(x3 - X2)0S) ix3, 62) + (a;2 + 3:3 - 2)0*(,) (x2, 62) 0i(a;3, ^2) 



■as(tj)/ij(xi, X2, X3, 61, 62) exp[-5'e/(tj)] 
- {a;20r,(a;3, &2)0^(') (a;2, ^2) 

+r^r^(,) (x2 + X3)0^(,) (x2, 62) + (3^2 - 2;3)0*(/) (a;2, ^2) 0^(a;3, ^2) 

+r^r^(,) (x2 - a;3)0j(,) (x2, ^2) + (a;2 + a;3)0*(/) {^2, 62) 0^(a;3, ^2) 
■a,(t))/i)(xi, X2, X3, 61, 62) exp[-S'e/(t))]} . 



(26) 



In the above equations, we have assumed that xi « X2,X3. Since the hght quark 
momentum fraction Xi in B meson is peaked at the small region, while quark momentum 
fraction X3 of 77 is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. 
After using this approximation, all the diagrams are functions of kf = xims/\/2 of B 
meson only, independent of the variable of . Therefore the integration of eq. fll3l) is 
performed safely. 

For the B — > rji]' decay, besides the Feynman diagrams as shown in Fig. 1 where 
the upper emitted meson is the r]', the Feynman diagrams obtained by exchanging the 
position of r] and rj' also contribute to this decay mode. The corresponding expressions 
of amplitudes for new diagrams will be similar with those as given in Eqs. (11611 26|) . since 
the 1] and rj' are all light pseudoscalar mesons and have the similar wave functions. The 
expressions of amplitudes for new diagrams can be obtained by the replacements 



For example, we find that: 



V 



(27) 



F , 



erji 



Farj' Far) i 



arj' 



PI 

ar] 1 



^P2 
arj' 



ar) ■ 



(28) 
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B. Mixing of t] and rj' meson 



Before we write down the complete decay amplitude for the studied decay modes, we 
firstly give a brief discussion about the rj — r]' mixing and the gluonic component of the rj' 
meson. There exist two popular mixing basis for rj — rj' system, the octet-singlet and the 
quark flavor basis, in literature. Here we use the SU{3)f octet-singlet basis with the two 
mixing angle scheme [2^, sO] instead of the simple one mixing angle scheme to describe 



the mixing of rj and r)' mesons, since the former scheme can archive a better agreement 
with the relevant data, such as the decay width of i]^'^ 77, the rj'^'^'j transition form 



factors, the radiative J/^ decays[30|, |3l[. In the two mixing angle scheme, the meson 77, 
1]' and the decay constants can be defined as 



r]\ ^ / cos 6*8 - sin 9i\ f m 
7]' J \ sin 6*8 cos 9i J 



(29) 



with the flavor S'f/(3)-octet and -singlet components 

?7i = {uu + dd + ss) I v^, r78 = {uu + rfd — 2ss) / -\/6, (30) 
in the quark model, and the two mixing angles Q\ and ^g, in principle, can be determined 



by various related experiments [30l . l31l| . In the numerical calculations, we will use the 
following 1] — 1]' mixing parameters 3^ 



^8 = -21.2°, 9i = -2.4°, 

/i = ISlMeV, /8 = 169MeV, (31) 



obtained by setting fg = f^ = 130 MeV and = v^S/l- - = 1.41/^ [SOj. The second 
set of the t] — rj' mixing parameters obtained phenomenological and used frequently in 
literature reads |i30'] 

^8 = -21.2°, 01 = -9.2°, 

/i = 1.17^ = 152MeV, /8 = 1.26^ = 164MeV, (32) 

for = ISOMeV. We usually use the first set of mixing parameters in numerical calcu- 
lations, unless otherwise stated. 

As shown in Eq. fl30l) . 1] and rj' are generally considered as a linear combination of light 
quark pairs uu, dd and ss. But it should be noted that the rj' meson may have a gluonic 
component in order to interpret the anomalously large branching ratios of B —>■ Krj' 
HQ m, 113. Although some progress have been achieved in recent years about this 



problem, we currently still do not know how to calculate reliably the gluonic contributions 
to the B meson two body decay modes involving rj' meson. For the studied decay modes in 
this paper, on the other hand, currently available measurements for the branching ratios 
still have big uncertainties. It is therefore reasonable at present to consider only the 
dominant contributions from the quark contents of f]' meson, while take the subdominant 
contribution from the possible gluonic content of r]' meson as a source of theoretical 
uncertainties. Following Ref. |33(], on the other hand, we also estimated the possible 
gluonic contributions to B ^ 'qi')fji') decays induced by the gluonic corrections to the B 
T]^'^ transition form factors [33] and found that these corrections to both the branching 
ratios and CP violating asymmetries are indeed very small. 
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C. Complete decay amplitudes 

For —>■ rji] decay, by combining the contributions from different diagrams, tlie total 
decay amplitude can be written as 



'7 5 2 1 1 1 1 

—^t { gC'a + -^Ci — 2C5 — -Cq — -C7 — -Cs + -Cg — -Cio 

—^t ( C-s + -C4 — C5 — gCe + + qCs — 2^*9 — g^io j 



+M,,F,{e,,es 



-6(C4--Cio)F2(^l,08j 



"6 ( C4 — -Cio ) -F2(^^i, ^s; 



1 



iuC2 — it[C^ + 2C4 — -Cg + 2^10 



2C6 + -Cg ) 9^) 



Cq - -Cg ) -^2(^1, 6*8)^ 



a?? 



6 



(33) 



where = V*fyud, = Vtl^ta, while the mixing parameters and the relevant decay 
constants are 



rd 



^cos^s-^^sin^^i, ^2(^1, 08 ) 



^sin08 + y^cos^i, (34) 



fi 



fi 



_ cos 9s — ^ sin 6'i, = — ^ cos ^8 + sin 9i. 
v6 v3 v3 v3 



/i 



(35) 



It should be mentioned that the Wilson coefficients Ci = Ci(t) in Eq. (1551) should be 
calculated at the appropriate scale t using equations as given in the Appendices of Ref . . 
Here the scale t in the Wilson coefficients should be taken as the same scale appeared in 
the expressions of decay amplitudes from Eqs. (ITB]) to (12^ . This is the way in pQCD 
approach to eliminate the scale dependence. 
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Similarly, the decay amplitude for 7777' can be written as 



(75 21111 
3C3 + ^C, - 2C, - -Ce - -Cr - + -C, - -Cio 

• I C'i -\- — Ca — C5 — ~C(\ ~\~ —C7 ~\~ ~Cr — —Co — — Cin 

st^dTg4 5 gb ^ i g« 2^ g^'J 
- ^s)/^^, + F^^^,Fi{9„ 9s) fi) 6 f + - ^Cy - 



+ (M,, + Me,0i^i(^i,^8)i^i'(^i,^8 



1 



[M,,F,i9^, 9s)F^i9u 9s) + M,^,F[{9^, 9s)F^{9^, 9s)] it[C^- -Cio 



- {M^^F,{9,, 9s)F^{9,, 9s) + M^^,F[{9,, 9s)F2{9,, 9s)) it [C, - 

+ (Mar, + Mar,') F^{9i, 9s)F[{9u 9 



iuC2 — it [ Ci + 2C4 — -Cg + 2^10 



- {Mar, + Mar,') F2{9^. 9s)F!,{9^, 9s)it C4 



1 



-C 
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- + M^, ) ^^8)F{(ei, 9s)it (c, - \ 



- (m5 + m5 ) F,{9u 9s)Fi{9u 9s)it ( 2C, + 

1 



- {M^^ + M^^, ) ^8)^2(^1, Os)it [Ce - -^8 

-/b • (Ff/ + 9s)F{{9u 9s)it + - - , (36) 

where the relevant mixing parameters and decay constants are 

Fii0i,08) = y^sin^s + y^cos^i, F^{9,,9s) = -^sm9s + ^cos9u (37) 

fi = sin ^8 + cos 9i, = — ^ sin 6*8 + cos 9i. (38) 
v6 v3 v3 v3 

The complete decay amplitude A4{r]'T]') for i? ^ t]'-)]' decay can be obtained easily 
from Eq. (!33|) by the following replacements 

/d rs ^ rd rs 
r, J r] J r]' J r]' 

Fi{9i,9s) F[{9^,9s), F2{9u9s) F^{9u9s). (39) 
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Note that the contributions from the possible gluonic component of r]' meson have not 
been included here. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate the branching ratios and CP violating asymmetries for 
those considered decay modes. The input parameters and the wave functions to be used 
are given in Appendix [Bl In numerical calculations, central values of input parameters 
will be used implicitly unless otherwise stated. 

Based on the definition of the form factor F^^^ as given in Eq. (fT7|l and consider 
the relation of = Fe^/ at the leading order of pQCD approach, we find the numerical 
values of the corresponding form factors at zero momentum transfer: 

F^riq' = 0) = ^^'(^^ = 0) = ^ = 0.30l°:°^K), (40) 



m 



B 



for Ub = 0.40 ± 0.04 GeV, which agrees well with those as given in Refs. [3J, [35 



A. Branching ratios 



Using the decay amplitudes obtained in last section, it is straightforward to calculate 
the branching ratios. For —>■ 7717, and rj'ri' decays, the decay amplitudes as given in 
Eqs. ( !33ll and (1361) can be rewritten as 



M = V:,V^,T - V*V,,P = V:,V^,T [1 + 



ze 



i{a+5) 



(41) 



where 







p 






T 



(42) 

is the ratio of penguin to tree contributions, a = arg —^^^^ is the weak phase (one of 

the three CKM angles), and 6 is the relative strong phase between tree (T) and penguin 
(P) diagrams ^. In the pQCD approach, the ratio z and the strong phase 6 can be 
calculated perturbatively. 

From Eq. ( |4T1) . it is easy to write the decay amplitude for the corresponding charge 
conjugated decay mode 



M = Vu,V:,T - Vt,V:,P = Vu,V:,T [l + ze^(-"+^)] . 
Therefore the CP-averaged branching ratio for B'^ rj^'^^'^ decays is 

Br = {\M\^+\M\^)/2= iKfeKd^l' [1 + 2^ cos a cos 5 + z^] 



(43) 



(44) 



^ The "T" and "P" term refers to the part of the decay amphtude Ai in Eq. (|¥T|) . which is proportional 
to V*f,Vud and V^lVtd respectively. 
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where the ratio z and the strong phase 5 have been defined in Eqs.f l41l) and ( l42i) . 

By employing the two mixing angle scheme of rj — t]' system and using the mixing 
parameters as given in Eq. (!3T]) . one finds the CP averaged branching ratios for the 
considered three decays as follows 

Br{ ^ T]T]) = [1.5 ± 0.3{uJb)to^2imo) ± 0.4(a)] x 10~^ (45) 
Br{ 5° ^ rjr]') = [omtoHitOb) ± 0.05(mJ) ± 0.04(a)] x 10"^ (46) 
Bri ^ r^V) = io.68+°:°^(a;,) ± 0.05(m^)l°:22(a)] x 10"^ (47) 

where the main errors are induced by the uncertainties of uj^ = 0.4 ± 0.04 GeV, rriQ = 
1.4 ± 0.1 GeV and a = 100° ± 20°, respectively. 

Besides the theoretical uncertainties as shown in Eqs. (145114 7p . the change of the mixing 
scheme oirj — rj' system or changes of mixing parameters within a given scheme may also 
induce moderate or even significant changes to the theoretical predictions. The central 
values of the branching ratios for the considered three decays, for example, will be 

Br{ B^ riri) \cv = 1-9 x 10"^ (48) 
Br{ B^ riri') \cv = 0.82 x 10"^ (49) 
Bri VV) \cv = 0.66 x 10"^ (50) 

if the second set of mixing parameters in the two mixing angle scheme, as given in Eq. (l32l) . 
are used, and 



Br{ B^ riri) \cv = 1-9 x 10"^ (51) 
Br{ 5° ^ riri') \cv = 1-2 x 10"^ (52) 
Br{ 5° ^ r/V) lev = 0.69 x 10"^ (53) 

if the one mixing angle scheme with 6p = —12.3° are employed. From the above numerical 
results one can see that (a) the induced variation is about 2% or 20% for B rj'rj' and 
B rjT] decay respectively, which is consistent with the general expectation; (b) for 
B 7]7]' decay, however, the resultant change of the branching ratio is roughly a factor 
of two. Such variation is channel-dependent and should be considered as one kind of 
theoretical uncertainties. The large change in Br{B rji]') is induced by destructive 
interference between the individual decay amplitudes from different Feynman diagrams 
when two mixing angle scheme are employed. 

As a comparison, furthermore, we also list here the theoretical predictions for the 
branching ratios of the three decays in the QCDF approach[5| and the SCET approach ^ 



^ The theoretical predictions for the branching ratio and the direct CP violating asymmetries in the 
SCET approach, as shown in Eqs. (|54ll56p and Eqs. (|70ll72|l . are directly quoted from Table VII of 



Ref. 



ITj . for the case of Theory II. 
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FIG. 2: The a dependence of the branching ratios (in units of 10 ®) 
(a) is for ml = 1.4 GeV, uj^, = 0.40 ± 0.04 GeV; and (b) is for Uf, = 
GeV. 



of rjrj decays. Here 

0.4 GeV, ml = 1.4 ± 0.1 
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Br{B^ ^ 7]r]) 
Br{B^ ^ r]r]') 
Br{ 5° ^ r/V) 



(0.16^°:^^) X 
(1.0 ± 1.5) X 

(0.16^°:?^) X 

(2.2 ±5.5) X 

(0.06^°:^?) X 
(1.2 ±3.8) X 



10 
10 

10 
10 

10- 
10- 



QCDF, 
SCET, 

QCDF, 
SCET, 

QCDF, 
SCET, 



(54) 
(55) 
(56) 



where the individual errors as given in Refs. and [13] have been added in quadrature. 
It is easy to see that (a) the theoretical predictions for Br{B — >■ rjr]) and Br{B rj'r]') 
in both the QCDF and the pQCD approaches agree very well; (b) for Br[B —>■ rjri'), the 
central value of the pQCD prediction is about half of the QCDF prediction, but still agree 
within one standard deviation; and (c) the central values of the theoretical predictions 
in the SCET approach are much larger than those in QCDF and pQCD approaches, but 
still consistent with them if one takes the very large theoretical uncertainties in SCET 
approach into account. 

It is worth of mentioning that the FSI effects are not considered here. The smallness 
of FSI effects for B meson decays into two light mesons has been put forward by Bjorken 
based on the color transparency argument 21], and also supported by further renor- 



malization group analysis of soft gluon exchanges among initial and final state mesons 

At present, the pQCD predictions still have large theoretical errors ( say ~ 50% ) 
induced by the large uncertainties of many input parameters. In our analysis, we con- 
sidered the constraints on these parameters from analysis of other well measured decay 
channels. For example, the constraint l.lGeV < ml < 1.9GeV was obtained from the 
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FIG. 3: The same as Fig. [2] but for rjrj' decay. 




FIG. 4: The same as Fig. [2] but for r\r\ decay. 



phenomenological studies for B ^ -kit decays |8|, while the constraint of a = (93lg^)° 



are obtained from the direct measurements [37|]. In estimating the uncertainties, we still 



conservatively consider the range of a = 100° ± 20°. 

In Figs. [2], [3] and HI we show the parameter dependence of the pQCD predictions 
for the branching ratios of B ^ rji] , rji]' and ri'rj' decays for Ub = 0.4 ± 0.04 GeV, 
rriQ = 1.4 ± 0.1 GeV and a = [0°, 180°]. From the numerical results, we observe that the 
pQCD predictions are sensitive to the variations of Ub, rriQ and a. 
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B. Branching ratios with updated wave functions 



In pQCD approach, the only input matters are wave functions of the involved particles, 
which stand for nonperturbative contributions. Since currently available wave functions 
are not exactly determined, large error of the theoretical predictions will be produced 
from the uncertainty of the relevant B, t] and rj' wave functions ^. 

In last subsection, we choose the pion DAs derived from QCD sum rule as given in 



Refs. |3J, |38| with a fixed decay constant = 190 MeV(see Appendix [B] for more 
details). The resultant pQCD predictions for CP averaged branching ratios are consistent 
with the measured values, which may be regarded as an indication that above inputs are 
reasonable. 

In this subsection, in order to check the theoretical uncertainty induced by the variation 
of wave functions, we recalculate the CP averaged branching ratios by employing the 
updated models of the pion DAs as given in Ref. [35|. Although the structure of 0^'^'^ 
as given in Appendix [B] remain unchanged, but the Gegenbauer moments and in 
the updated pion DAs are changed significantly, they are now 

{al,al) = (0.115,-0.015), (57) 

instead of the old (0.44, 0.25) as shown in Eq. (]B5p . 

Using the updated pion DAs with (a^,aj) = (0.115, -0.015) and fs = 210 MeV |39 
we find numerically that 

(i) The values of the form factors as given in Eq. (I40i) will be decreased by about 
10%, and we now have 

Fo'^ri^" = 0) = Foy^'iq' = 0) = 0.27l°:°^(c.,), (58) 
for Ub = 0.40 ± 0.04 GeV. 

(ii) The CP averaged branching ratios of the considered decays are 

Bri B' vv) = [0.67i°:?^(a;,)l°:°^(mJ);°i^(a)l°:J^(4)] x 10"^ (59) 
BriB'^W) = [0.18;°:°^(o;,)±0.01(mJ)±0.01(«)l[J:J°(4)] xlO-^ (60) 
BriB'^vW) = [0.1i;°:°t(c^,)±0.01(mJ)±0.04(a)l°:J°(a^^^ (61) 

where the main errors are induced by the uncertainties of Ub = 0.4 ± 0.04 GeV, 
= 1.4 ± 0.1 GeV, a = 100° ± 20°, = 0.115 ± 0.115, respectively. The errors 
induced by varying aj in the range of a1 = —0.015 ± 0.015 are very small: only 
about ±0.03, ±0.004 and ±0.01 (in unit of 10"'') for B rjrj^rjr]' and rj't]' decay, 
respectively. We here have assumed that the updated Gegenbauer moments and 
a1 can vary by 100%. 



3 



This constitutes a large theoretical uncertainty common to the pQCD, QCDF and other factorization 
approaches where the meson wave functions are used as input. 
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FIG. 5: and aj-dependence of the CP averaged branching ratio (in units of 10 Br{B 
rj-q) for < < 0.50 (a) and -0.10 < aj < 0.30 (b). 

It is easy to see that the pQCD predictions in Eqs. fl59ti6T|) is smaller than those in 
Eqs.f H5M7|) by roughly a factor of 2 — 5. It is not difficult to understand such situation: 
the large decrease is induced by the great changes of the updated Gegenbauer moments 
in the wave functions. The updated Gegenbauer moment 03 = 0.115 is indeed much 
smaller than the previous one 0.44 for the leading twist-2 distribution amplitude, a factor 
of four decrease in magnitude. For Gegenbauer moment aj, which governs the high order 
contributions in pion DA, it changes from = 0.25 to = —0.015, even the sign is 

changed yet. The form factor Fq^^ then reduces from 0.30 to 0.27, which leads to a 

smaller branching ratio. Besides the effect of a smaller form factor F^[^^^ ^ which were 
extracted from the Fig. 1(a) and 1(b), the contributions to the branching ratios from 
the remaining Feynman diagrams l(c)-l(h) also become smaller than before due to the 
significant variations of 03 and 04. The total decrease of the branching ratios therefore 
becomes significant. This fact tell us that current theoretical predictions still have a 
strong dependence on the form of meson wave functions. 

We now take B ^ rjr] decay as an example to show the Og- and -dependence of the 
theoretical predictions explicitly. Firstly, if we vary only the value of Gegenbauer moment 

in the range of < ^ 0-50, the pQCD prediction will be 

0.55 X 10"^ < Br{B r]r]) < 1.23 x 10"^. (62) 

This can be seen more clearly in Fig. [5^, where the CP averaged branching ratio Br{B 
Tjr]) shows a linear dependence on Og. For the parameter aj, we also find the similar linear 
dependence: 

0.53 X 10"^ < Br{B r]r]) < 1.40 x 10"^ (63) 
for -0.10 <al< 0.30. 
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Although the pQCD predictions for the branching ratios obtained by employing the 
previous or updated pion DAs are all consistent with the measured values due to still 
large theoretical and experimental errors, the pQCD predictions in Eqs. fl59tl6T|) are more 



are 



reliable, in our opinion, since the updated Gegenbauer moments as given in Ref. [35 
used to obtain these results. Of course, better wave functions of particles involved are 
clearly needed to reduce the errors of theoretical predictions. 

In next section, we will use the updated Gengenbauer moments {a^, al) = 
(0.115, —0.015) to calculate the CP-violating asymmetries for the three considered de- 
cay modes. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of B ^ decays 
in pQCD approach. We here use the wave functions as presented in Appendix [Bl but with 
the updated Gengenbauer moments (02,04) = (0.115, —0.015). Because these decays are 
neutral B meson decays, so we should consider the effects of B^ — B^ mixing. For B^ 
meson decays into a CP eigenstate /, the time-dependent CP- violating asymmetry can 
be defined as 

Br (B°{t) ^ f) - Br {B^{t) ^ /) 

^ ^ - A'^SpCos{Amt)+ASpsm{Amt), (64) 



Br (B\t)^ f)+Br{B^{t)^ f) 



where Am is the mass difference between the two B^ mass eigenstates, t = tcp — Uag 

is the time difference between the tagged B^ (B ) and the accompanying B (5°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
mixing induced CP-violating asymmetries A^p and A^p can be written as 



- i+iAcpr '''' ~ i+iAcpr ^ ^ 

where the CP- violating parameter \cp is 

Here the ratio z and the strong phase 5 have been defined previously. In pQCD approach, 
since both z and b are calculable, it is easy to find the numerical values of A^^p and ACJ^p 
for the considered decay processes. 

By using the mixing parameters in Eq. (l3T|l and the input parameters as given in 
Appendix B, one found the pQCD predictions (in units of 10~^) for the direct and mixing 
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induced CP-violating asymmetries of the considered decays 

^CP - VV) = +53.5tn(«)ll7(^.)loi«), (67) 

AdFiB' - W) = -13.lli^i(«)i^:^(c.,)™(rnJ), (68) 

AfpiB'^r^W) = +23.7™(«)iil^(^6)l|^K), 

^—(5° ^ r/V) = +93.2^|9(a)^?i2i(a;,)l|?K), (69) 

where the dominant errors come from the variations of a = 100° ± 20°, Ub = 0.4 ± 0.04 
GeV and = 1.4 ± 0.1 GeV. 

In Fig. El we show the a— dependence of the pQCD predictions for the direct and the 
mixing- induced CP-violating asymmetry for — ^ rjr] (dotted curve), B^ r]r]' (solid 
curve) and B^ — > rj'r]' (dashed curve) decay, respectively. 

As a comparison, we present the theoretical predictions for Aqp^B'^ — >■ r]^'^r]^'^) (in 
units of 10~^ ) in both the QCDF approach [5(] ^ and the SCET approach [13] 



Adir (T^O \ _ f +63^74, QCDF, . . 

■^cp[^ ^VV) - I +48 ±32, SCET, ^^^^ 

/idiri'nO _ / +561^44, QCDF, . ^ 

^cpl^ ^m) - I +70 ±24, SCET, ^'^^ 
■Acp[^ ^VV) - I +60 ±38, SCET, ^^^^ 

where the individual errors as given in Refs. jH| and [l3] have been added in quadrature. 
From above numerical results one can see that 

(i) In the considered three kinds of factorization approaches, the theoretical pre- 
dictions for the CP violating asymmetries are generally large in magnitude, and 
consistent with each other if one takes the very large theoretical uncertainty into 
account. 

(ii) For B rji] decay, the pQCD prediction for A^p has an opposite sign with 
those given in QCDF and SCET approach, which may be tested in the future B 
experiments. 

If we integrate the time variable t, we will get the total CP asymmetry for B'^ ^qid^qi') 
decays. 



There is a sign difference between the term A^^p defined here and the term Cf defined in Ref. i.e., 

Adir 
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FIG. 6: The direct and mixing-induced CP asymmetry (in percentage) of — > rj-q (dotted 
curve), r]!]' (solid curve) and r/' 77' (dashed curve) decay as a function of CKM angle a. 



where x = Am/T = 0.771 for the — B mixing Numerically, we found (in units 
of 10-2) that 



A'Sp{B^ 
A'SpiB"^ - 



rjT]) 
> r]r]') 

v'v') 



42.2l^^:?( 



v+4.0/ , N+3.2, 



m. 



59.9l°:°(«)l|g(c.,)^i«(m2 



(74) 
(75) 
(76) 



D. Effects of possible gluonic component of r]' 



Up to now, we have not considered the possible contributions to the branching ratios 
and CP- violating asymmetries of S — r]^'^r]^'^ decays induced by the possible gluonic 
component of 77' (isl . [3, 41|- When a non-zero gluonic component exist in rj' meson, an 
additional decay amplitude Ai' will be produced. Such decay amplitude may interfere 
constructively or destructively with the ones from the qq {q = u, d, s) components of rj', 
the branching ratios of the decays in question may be increased or decreased accordingly. 

In Ref. |32[], Beneke and Neubert computed the leading two-gluon contribution to the 
B 7]^'^ form factors using the framework of QCD factorization. In Ref. T^, Kou 
examined the gluonic component of rj' and the contributions to the process gg — >■ t]'. In 



his paper [13|, the rj and rj' meson were written as 



\r]> = Xr,\r]g > +Yri\r]s >, 
\ri' > = Xri'\rig > +Yr)i\ris > +Zn'\gluonium >, 



(77) 



where r]g = {uu + dd) / \/2 and rfg = ss. From the experimental data on the radiative 
light meson decays, such as ^ v'l^v' ^ {4>^P:l)l J/\psi — > ri''y decays, the author 
found that the gluonic component in r]' should be less than 26%. 

In Ref. by employing the pQCD factorization approach, Charng, Kurimoto and Li 
calculated the flavor- singlet contribution to the B — ^ r]^'^ transition form factors induced 
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by the Feynman diagrams with the two gluons emitted from the hght quark of the B 
meson (see Fig. 1 of Ref. [s^]), and they found that this gluonic contribution is neghgible 
(~ 5%) in the B ^ rj form factors, but can reach 10% — 40% in the B ^ rj' ones in the 
whole parameter space. Such enhancement to B ^ rj' transition form factor can be help 
to explain the large branching ratio Br{B Krj'). 

In order to check the gluonic effects on the decay modes under study, we here follow 
the same procedure as being used in Ref. [s^] to include the possible gluonic contributions 
to the B 7]^'^ transition form factors Fq^^ and in turn to the branching ratios and 
CP violating asymmetries. We found that the gluonic contributions to i? ^ r]^'\^'^ decays 
are rather small: 

• The central values of the pQCD predictions for the branching ratios remain basically 
unchanged ( the variation is less than 2% ) for B rji] and B Vjvl decay. For 
B r]r]' decay, the enhancement is only about 5%: the central value is changed 
from 0.60 x 10"^ to 0.63 x 10"^ 

• As for the CP-violating asymmetries of B ^ 7^(0^(0 decays, the possible gluonic 
corrections are largely canceled in the ratio, and therefore negligible (less than 3%). 

• The smallness of the gluonic corrections to the branching ratios can be understood 
as follows: (a) the gluonic correction to B ^ rj transition form factor itself is 
negligibly small [3^; (b) only the first two diagram Fig. 1(a) and 1(b) are affected 
by the gluonic corrections to -B — > 77' form factor, while the contributions from other 
six diagrams remain unchanged; and (c) the total effects are consequently small. 

Although much progress have been achieved in recent years, but frankly speaking, 
we currently still do not know how to calculate reliably the contribution of the possible 
gluonic component in rj' meson. From our previous works, as presented in Refs. [ll,|l3 



where only the dominant contributions from quark contents of r] and 77' were taken into 
account, the pQCD predictions for the branching ratios of -B — > pr]^'^ and B — >■ tit]^'^ 
decays also show a very good agreement with currently available data. It seems that large 
gluonic contributions are unnecessary for these decay modes. For B Krj' decays, on the 
contrary, th e g luonic contribution may play an important role in explaining the so-called 



Kf]' puzzle [32, l33 



Of course, more theoretical studies about the effects of possible gluonic component 
in r]' and better experimental measurements for the relevant decay modes are needed to 
clarify this point. 



V. SUMMARY 

In this paper, we calculated the branching ratios and CP-violating asymmetries of 
B^ — »• 7777, 7777' and r^'rj' decays at the leading order by using the pQCD factorization ap- 
proach. Besides the usual factorizable diagrams, the non-factorizable and annihilation 
diagrams are also calculated analytically in the pQCD approach. Furthermore, the anni- 
hilation diagrams provide the necessary strong phase required by a non-zero CP- violating 
asymmetry for the considered decays. 

From our calculations and phenomenological analysis, we found the following results: 
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The pQCD predictions for the form factors of B — » ?7 and t]' transitions agree well 



with those obtained in QCD sum rule calculations [3J, |35 



Using the two mixing angle scheme and the updated Gegenbauer moments 03 and 
, the pQCD predictions for the CP-averaged branching ratios are 

Sr(i?° r]r]) = (0.67l°:i) x 10"^ (78) 
Br{B°^riri') = (0.18 ± 0.11) x 10"^ (79) 
Br{B'^v'v) = (0.1ll°;S^) X 10-^ (80) 

where the various errors as given in Eqs. (159116 ip have been added in quadrature. 
The leading pQCD predictions are consistent with currently available data, but 
both the theoretical and experimental errors are still large. 

For the CP-violating asymmetries of the considered three decay modes, the pQCD 
predictions are generally large in magnitude, and have large theoretical uncertainty. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function /ij's, coming from the Fourier transformations of the function 

he{xi,X3,bi,b3) = Kq {^Jxix^niBbi) [9{bi - b3)Ko {^/xiniBbi) Iq {^/xlmsb^) 

+e{bs - bi)Ko (v^msfeg) Jo {^M^Bbl)] St{x3), (Al) 

haix2, X3, 62, ^3) = Kq {iy/x2X3mBb2) [9{bs - 63)^0 (^v^^^B^s) Iq {iy/^ruBbi) 

+e{b2 - bs)Ko (^v^mfifea) lo (^v^msfog)] Stix^), (A2) 

hf{xi,X2,X3,bi,b2) = 1 61(62 - bi)lo{MBy/XiX3bi)Ko{MBy/XiX3b2) 

>| / Ko(MbF(i)62), for Ffi)>0\ 
+ (&i ^b2)j- |^^H«(M,^ 62), for F^,, < oj ' ^^3) 

/iJ(xi,X2,X3,6i,62) = ^0{bi - b2)Ko{iy/x2X3biMB)lo{iy/x2X3b2MB) + {pi ^ 62) I 

■Kq^^/xi + X2 + X3- X1X3 - X2X3 6iMb), (A4) 
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hj{xi,X2,X3,bi,b2) = < 6'(6i - 62)Ko(iVxia^6iMij)Io(zyx^fe2MB) 



+ {bi 



Ko(MbF(2)6i), for F2)>0 



f H«(Ms I 60, for <0 



, (A5) 



(2) 



(2) 



where Jq is the Bessel function and Kq, Iq are modified Bessel functions Kq^—ix) 
— {-n /2)Yq{x) + z(7r/2) Jo(a;), and i^(j)'s are defined by 



1^ 

^(2) 



(Xi - X2)X3 . 



The threshold resummation form factor St{xi) is adopted from Ref . |28 

2i+2T(3/2 + c 



5"/ fx) 



v^r(i + c) 



[x(l -x)]' 



where the parameter c = 0.3. This function is normahzed to unity. 
The Sudakov factors used in the text are defined as 



Sabit) 



Scd{t) 



s [xiniB/ 6i j + s [x^niB/ v^, 63) + (^(1 - x-^wlbI V2, k 
1 



-ln(6iA) -ln(63A) 
s {xiitlb/ V2, 61 j + s {x2mB/ V2, b-^ + s ^(1 - X2)mB/ v^, h 

s I (1 - x-i)niB/V2,bi 
ln(t/A) 



+s {x^itlb/V^, bi 
1 



, ln(t/A) 
2 In + In ■ 



Sefit) 



ln(6iA) ' '"-ln(62A) 
s {xiTUb/ V2, 61 j + s {x2mB/ V2, 62 j + s (^(1 - X2)mB/ V2, &i 

+s ^X3mB/ v/2, 62) + s (^(1 - x^)mB/V2, 62 j 

ln(t/A) 



1 

"A 



, ln(t/A) 
In + 2 In ■ 



Sgh{t) 



ln(6iA) -hi(62A) 
= s {x2mB/ v^, 62) + s [x^niB/ v^, 63) + s [il 

+ s (^(1 - xz)mB/ V2, ^3) - ^ 



, ln(t/A) , 
In . /. .\ + In 



X2)mB/V2, 62 
ln(t/A) 



(A6) 
(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



(A12) 



-ln(62A) -ln(63A)J ' 
where the function b) are defined in the Appendix A of Ref.jl]. The scale tj's in the 
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above equations are chosen as 

tl = max{y/x^mB, l/h, I/63) , 
tl = max(yxrmB, 1/bi, I/63) , 
tl = max{y/x^mB, I/62, l/h) , 
tl = max{^/x^mB, 1/62, 1/&3) , 
tf = maxi^y/xiXsniB, a/ (xi - X2)x^mB, I/&1, 1/&2) , 

= max(Va;i + X2 + X3 - Xia;3 - X2X3mB, V^2a;3'^B, l/^i, I/&2) , 
tj = max{y/x^mB, - X2)x3mB, l/h, I/62) • (A13) 

They are chosen as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 



APPENDIX B: INPUT PARAMETERS AND WAVE FUNCTIONS 



In this Appendix we show the input parameters and the light meson wave functions to 
be used in the numerical calculations. 

The masses, decay constants, QCD scale and meson lifetime are 



A. 



MS 



250MeV, U = ISOMeV, /b 



m, 



1.4GeV 



m, 



riss 



5 



2.4GeV, 



190MeV, 

= 140MeV, fx = leOMeV, 



1-12^ 



Mb = 5.2792GeV, Mw = 80.41GeV, tbo = 1.54 x 10" "s 
The central values of the CKM matrix elements as given in Ref. lii are 



\Vud\ 

\Vtb\ 



0.9745, IK 



ub I 



0.0040, 



0.9990, \Vtd\ = 0.0075. 
For the B meson wave function, we adopt the model 



(f)Bix,b) = Nbx'^{1 — xYexp 



Ml x" 



Mb)' 



(Bl) 



fB2l 



(B3) 



where Uf, is a free parameter and we take Uf, = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for Ub = 0.4. This is the same wave functions 
as being used in Refs. [8|, iS, |28i] , which is a best fit for most of the measured hadronic B 
decays. 



For the distribution amplitudes 



Vdd' ^Vdd 



_ and appeared in Eq. (ITll) . we utilize the 
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result from the light-cone sum rule 34, 13] including twist-3 contribution: 



-^/..a;(l - a;) <j 1 + a^^ [5(1 - 2xY - l] 



15 



[21(1 - 2xY - 14(1 - 2xf + 1] I 



X) 



1 
8 
3 



/. 1 + 



+ 



I (3°"' - f "L) [3(1 - ^-r- - 1] 

[35(1 - 2x)^ - 30(1 - 2xy + 3] 



2 
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V2K 
1 

6^ 



/x(l - 2x) 



{5V3 - - ^p',, - \p\/f)i\^x^ - lOx + 1) 



with 34, 38 



J')dd 



PVd 



al = 0.44, 



-.Vdd 



0.25, 



= 0.015, Us = -3.0. 



The updated Gegenbauer moments as given in Ref. [35 



are 



0.115, 



-0.015. 



(B4) 



(B5) 



(B6) 



We also assume that the wave function of uu is the same as the wave function of dd 



14l |. For the wave function of the ss components, we also use the same form as dd but 
with uiq'^ and L instead of mf^ and fx, respectively. For and fy, we use the values as 
given in Ref. [38| where isospin symmetry is assumed for and SU{3) breaking effect is 
included for fy-. 



= fy=^Vl-fl (B7) 

These values are translated to the values in the two mixing angle scheme, which is 
often used in vacuum saturation approach as: 

/i = 151MeV, /s = 169MeV, 



-2.4°, e. 



-21.2°. 



The parameters uiq {i = '/^^^(u^i), r/ss) are defined as: 



m, 



Vdd(uu) 





{rriu + ma) 



2Ml 



{2ms 



(B8) 



(B9) 
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